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Abstract. The technique of adjoint cascade equations has been applied to calculate 
the properties of air showers produced by extremely high energy (EHE) 7-rays in the 
energy range f 18 — f 22 eV. The high intrinsic accuracy, combined with very modest 
(compared with the traditional Monte Carlo codes) computational time requirements, 
make this method as an effective tool for detailed study of development of EHE showers 
in the Earth's atmosphere. In this paper a wide range of parameters of 7-ray induced 
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l/") . showers are analysed taking into account two independent effects which become crucial 

for the cascade development in the EHE regime - the Landau-Pomeranchuk-Migdal 
(LPM) effect and the interaction of primary 7-rays with the geomagnetic field (GMF). 
Although the LPM effect leads to dramatic modifications of shower characteristics, 
C"| . especially at primary energies exceeding 10 19 eV, the GMF effect, which starts to 

♦~h| "work" approximately at same energies, prevents, to a large extent, the LPM effect 

O ' by converting the primary 7-ray into a bunch of synchrotron 7-ray photons with 

-*— > ■ energies effectively below the threshold of the LPM effect. This bunch of the secondary 

C$ • photons hits the atmosphere and creates a large number of simultaneous showers. 

The superposition of these independent showers mimics a single shower with energy 

E = ££j ~ Eq, but without the signatures of the LPM effect. This makes the 

jJJ ■ longitudinal profile of such a composite electromagnetic shower quite similar to the 
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longitudinal profile of hadronic showers. At the same time, the number of muons as well 
as their lateral distribution differ significantly from the corresponding parameters of 
proton-induced showers. In the "gamma-ray bunch" regime, the total number of muons 
is less, by a factor of 5 to 10, than the number of muons in hadronic showers. Also, 
compared with the hadronic showers, the electromagnetic showers are characterized by 
a significantly narrower lateral distribution of muons. Even so, for inclined EHE 7-ray 
showers the density of muon flux at large distances from the shower core (> 1000 m) 
can exceed the electron density. 
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1. Introduction 

The question of the origin of most energetic particles observed in cosmic rays at 
energies E > 10 18 eV, often called Extremely High Energy cosmic rays (EHECRs), is a 
subject of intensive astrophysical speculations and debates. Within the conventional 
acceleration (or "bottom- up") scenarios, the powerful radiogalaxies, active galactic 
nuclei and clusters of galaxies (see e.g. |j], |2|, f|), as well as transient events like the 
sources of 7-ray bursts f|, f| || are believed to be the most probable sites of production 
of EHE cosmic rays. These models have, however, a difficulty connected with the 
large, typically > 100 Mpc distances to the most prominent representatives of these 
source populations. Indeed, if the observed EHE cosmic rays are protons produced 
in distant extragalactic sources, their interaction with the 2.7 K cosmic microwave 
background radiation (MBR) should give rise to the so-called Greisen-Zatsepin-Kuzmin 
(GZK) spectral cutoff. Depending on the characteristic distance scale to the ensemble 
of sources responsible for the bulk of the observed EHE cosmic ray flux, the position 
of the cutoff is expected between £~5x 10 19 and 10 20 eV (see e.g. Refs.|F], ||, |[). 
Therefore, the registration of the EHECR events close to 10 20 eV by Haverah Park, 
Fly's Eye, Yakutsk, and AG AS A detectors (for review see Refs.|I| |TT|) pose a serious 
challenge for any of these conventional models. This difficulty initiated a quite different 
approach to the solution of the problem of origin of EHE cosmic rays - the hypothesis 
of non- acceleration (or "top-down") scenario which assumes that the highest energy 
cosmic rays are result of decay of hypothetical massive relic particles originating from 
early cosmological epochs (see e.g. Ref.|12|] and references therein). A distinct feature 



of the "top-down" models of EHECRs is an unusually high content of 7-rays with the 
gamma/proton ratio > 1 at E > 10 20 eV [|l^, [14], 15|. Thus the gamma/proton ratio 



can serve as a conclusive diagnostic tool for observational proof (or rejection) of the 
"top-down" scenarios. 

Actually, a non-negligible content of 7-rays into the EHE cosmic rays is expected 
also in the conventional particle acceleration models. Within these models the highest 
energy 7-rays are contributed through the decay of secondary 7r°-mesons produced at 



interactions of EHECRs with 2.7 K MBR [[T(| 17|. If the average magnetic field in 
the intergalactic medium is less than 10~ 12 G, the development of the ultrarelativistic 
electron-photon cascades in the field of 2.7 K MBR makes the attenuation length of 
EHE 7-rays at E > 10 20 eV larger than the attenuation length of protons. Therefore, 
the 7-ray content in EHE cosmic ray flux depends strongly on the typical distances 
to the cosmic ray sources, and on the strength of the intergalactic magnetic field. It 
could be as large as 10 per cent in the case of homogeneous distribution of the EHECR 
sources in the Universe, but less than 1 per cent if the bulk of the observed EHECR 
flux is contributed by relatively nearby sources concentrated in our Local Supercluster 



of Galaxies. Thus the measurements of the gamma/proton ratio in EHE cosmic rays by 
forthcoming powerful instruments like the Auger Observatory or the Telescope Array 
(see e.g. Ref . fT0|, |TT|| ) would not only allow crucial tests to distinguish between the "top- 
down" and "bottom- up" scenarios, but also can provide important constraints on the 
spatial distribution of the EHECR accelerators within the conventional "bottom- up" 
models. 

The effective separation of 7-ray induced showers from the proton-induced showers 
is possible because of noticeable differences in the development of electromagnetic and 
hadronic cascades. Starting from energies of about 10 TeV, the low-energy muon content 



becomes a reliable criterion for recognition of electromagnetic showers [[18|, [19|, |2(| |21|). 
At energies E > 10 16 eV the content of muons in the electromagnetic showers increases 



with an increase of energy significantly faster than in the hadronic cascades |[22|| . The 
results of this paper based on detailed numerical calculations of characteristics of 
electromagnetic showers, using the so-called adjoint equation technique (see for review 
||23||), generally confirm this effect. However, the present detailed numerical study shows 
that the prediction of Ref. p2| concerning the absolute muon number, based on an 



approximate analytical approach, was significantly overestimated. 

At primary energies above 10 19 eV the Landau- Pomeranchuk-Migdal (LPM) effect 
|24| , |25|j leads to significant suppressions of the Bethe-Heitler cross-sections for the 
pair-production and bremsstrahlung processes in the Earth's atmosphere, and thus 
dramatically changes the character of development of electromagnetic showers. In 
practice, however, the interesting features of electromagnetic showers caused by the 
LPM effect do not appear in a "pure" form, but are largely compensated by another 
effect connected with interactions of EHE 7-rays with the geomagnetic field (GMF) 
26| . Remarkably, for the combination of relevant parameters, namely the strength and 



orientation of the geomagnetic field from one side, and the height, density and the 
mean atomic number of the Earth's atmosphere from another side, these two effects 
start to "operate" almost simultaneously, at energies above several times 10 19 eV. The 
interaction of 7-rays with the GMF by e + , e~ pair production does not imply, however, 
absorption of primary photons. In fact, the secondary (pair produced) electrons and 
positrons in the same geomagnetic field quickly produce a bunch of synchrotron photons 
with energies typically between 10 15 and 10 19 eV. The secondary photons initiate a large 
number of simultaneous electromagnetic cascades in the atmosphere. Because almost the 
whole energy of the primary photon is re-distributed between the synchrotron photons, 
the superposition of these showers mimics a single electromagnetic shower with total 
energy rather close to the energy of the primary particle. Moreover, because of the 
relatively low energies of the bunch photons, the LPM effect is essentially "switched 
off", and the longitudinal profile of the shower returns to the "standard" shape of 
electromagnetic showers without the signature of the LPM effect WM. For typical 



strengths of the GMF of about 0.2—0.6 G and for energies of primary 7-rays E > 10 20 eV, 
this interesting interference of two strong effects becomes unavoidable almost at all 
arrival directions of primary 7-rays. Therefore, any study of electromagnetic showers of 
such high energies requires thorough treatment of both the LPM and GMF effects. The 
depth profiles of the electronic component of EHE electromagnetic showers including 
these two effects have been studied by several authors |22], ^7], ^, |29|, [30]]. Recently the 



characteristics of the muon component have been calculated by Capdevielle et.al. [[H 
using the CORSIKA Monte Carlo code. They took into account the LPM effect, but 
ignored the interactions of the primary 7-rays with the geomagnetic field. We generally 
confirm the results of these studies, but significantly extend the range of calculated 
characteristics. In particular, we calculated such basic parameters of electromagnetic 
showers as the electron size, the depth profiles of cascade particles, the content and 
the lateral distribution of muons and compared them with the relevant parameters of 
hadronic showers. 

2. Interaction processes 

In calculations of characteristics of electromagnetic air showers we take into account 
the following interaction processes - the bremsstrahlung and the ionization losses for 
electrons and positrons; the pair production and the Compton scattering for photons. 



The LPM effect |2|], [25] is taken into account for the bremsstrahlung and the pair 
production processes. In addition, we take into consideration the interaction of EHE 
7-rays with the geomagnetic field. The cross-sections of magnetic cascade processes 
(the pair production and the synchrotron radiation) are calculated in accordance with 

We calculate the spatial distribution of the GMF within the framework of the dipolar 
central model (DCM; see e.g. |34)). In table 1 we present the GMF for two proposed 



sites of the Auger Observatory in Northern and Southern Hemispheres [35 1. Besides the 



DCM results we present in table 1 the strengths of GMF obtained in Ref.|36| by use of 



the data of the International Geomagnetic Reference Field |p7 |. One can see that the 
dipolar central model overestimates significantly the -B-field at the Southern Hemisphere 
site. Therefore in calculations we multiply the DCM results on B by special correction 



factors providing the coincidence with data of |36| near the Earth surface. With such a 



correction the field strength for the Northern Hemisphere site of the Auger Observatory 
is larger by a factor of ~ 2 compared with the Southern Hemisphere site (table 1). 

In figure 1 we show examples of the spatial profile of the GMF component (B n ) 
perpendicular to the arrival direction of the air shower (note that this component of the 
field is the only parameter which is needed to calculate the cross-sections of both the 
pair-production and synchrotron processes). The parameter t in figure 1 is the distance 



Table 1. The strength of the geomagnetic field (in Gauss) near the Earth surface 



Site location 


El Nuhuil, 


Millard County, 




Mendoza (Argentina) 


Utah (USA) 


Geographical coordinates 


35.2° S, 69.2° W 


39.1° N, 112.6° W 


Observation level, g/cm 2 


890 


870 


DCM data 


0.376 


0.497 


Data of @ 


0.250 


0.528 



between the Earth surface and the observation point measured along the air shower 
arrival direction [|. One can see that B n increases with the zenith angle 9 and decreases 
rapidly as a function of distance t. 

For calculations of basic characteristics of photonuclear interactions of 7-rays we 
performed special simulations with use of the SOPHIA generator p?J. This generator 
provides a comprehensive tool for simulations of photohadronic processes. In particular, 
this code includes the excitation and decay of baryon resonances, the direct production of 
pions and the multiple production of hadrons. To simplify the calculation technique we 
assumed that all unstable particles generated in the photoproduction process (except 
for muons) immediately decay at the point of their generation, i.e. we neglect the 
finiteness of life time of such particles ||. In this approximation all secondary muons 
are assumed to be emitted at the point of the primary 7-ray interaction. Thus, the 
inclusive differential (over the muon kinetic energy (T) and the muon emission angle 
($)) cross-section w il {E 11 T 1 'd) becomes the only input parameter. This cross-section is 
normalized as: 

2tt / oh? sin $ / w lJL {E 11 T ) d)dT = cr p h • n M (1) 

where a p h is the total cross-section of photoproduction, n M - the mean multiplicity 
of muons. Following to recommendations of Jl8| we define the quantity cr p h as 
Oph = o" 7P • A ' 91 where <t 7P is the total cross-section of 7-ray - proton interaction (we 
use for this quantity data from f38fl), A is the mean atomic number of air nuclei. 

Besides the electromagnetic cascades we have calculated also the characteristics of 
proton-induced showers. A standard quark-gluon-string (QGS) model is applied in this 



case to describe the hadronic interactions |33, 40 . 



f Note that for inclined showers the B-ficld depends also on the azimuth angle. However, for simplicity, 
for such showers we assume a fixed azimuth angle corresponding to the arrival direction from the nearest 
geomagnetic pole. 

f Our calculations show that this is a quite acceptable assumption, if one deals with muon detectors 
with sufficiently small threshold energy, < 1 GcV. 
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3. Technique of adjoint cascade equations 

The calculations presented below are obtained with the help of numerical solution of 
the adjoint cascade equations. In this section we describe shortly the main features of 
these equations|]. 

Let us consider as an example the system of adjoint equations that describe the 
longitudinal development of the electromagnetic cascade. Studying the longitudinal 
development of such a cascade one can neglect the lateral displacement of secondary 
particles and their angular distribution induced mainly by the multiple Coulomb 
scattering. The equation system has in this case the following form 

rE w ee (E,E')f(t,E')dE'-[ E 

lEo JE 



X M + f~ I" W ^ E > E')f(t, E')dE'- l* w ei (E, E')g(t, E')dE' = X e F,{2) 
at JE Q Jeo 

X i%+9-J^w ie {E,E')f{t,E')dE'~j\ 11 {E,E')g{t,E')dE' = A 7 G.(3) 



The adjoint functions / and g in Eqs. (|j), and @ describe the contributions 
(averaged over random realizations of the cascade development) to the response of some 
detector from a cascade generated by a primary electron (/) or a photon (g) of energy 
E; t is the distance between the observation level and a point of the primary particle 
appearance (measured along the particle arrival direction); E is the particle detector 
threshold energy; X a is the mean free pathlength for the cascade particle of type a 
(a = e or 7); functions w a p{E,E') define the differential (over the energy E') spectra 
of secondary particles of type (3 ((3 — e or 7). 

Properties of the detector are defined in Eqs. (|2|) and (Q) by the right hand side 
functions F and G and the boundary conditions f(t = 0,E) and g(t = 0,E). For 
example, if the detector measures at the observation level the total number of cascade 
electrons with energy E > Eq, then 

F(t,E) = G(t,E)=0; f(t = 0, E) = H(E - E ), g(t = 0,E) = (4) 

where H(x) = 1 for x > and H(x) = for x < 0. 

In the case of calculation of characteristics of muons produced in the electromagnetic 
cascade due to nuclear interactions of cascade photons one can use a system of adjoint 
equations similar to Eqs. (0) and ([3D, but having the different form of right hand 
side terms and boundary conditions. For example, if the detector registers the muons 
transversing a circle of a given radius (Rq) with the centre at the shower axis and 
orientated perpendicularly to this axis, we have 

f(t = 0,E) = g(t = 0,E) = 0, F(£,£) = 0, (5) 

I More detailed description of the adjoint equation technique is presented in |23||. 



G(t, E) = 2vrn J dd sin d J w,(E, T, #)p d (T, d)p m (T, d)H{f-T th )dT ,(6) 

where w^ is the differential cross-section of muon photoproduction defined in section 2, 
T t h = E — m^c 2 - the detector threshold kinetic energy, n - the concentration of air 
nuclei. The function 

p d (T) = exp{- costf" 1 / T [A dec (T') • ^T'^dT'} (7) 

Jt 

defines the probability for the muon to survive on the way from its production point to 

the observation level. T is the energy of the muon at the observation level found from 

the following equation 

t/ cos ■& = f dT'/e{T') , (8) 

e is the ionization loss rate, and Ad ec is the decay mean free pathlength. 

The function p m is the probability of the muon to escape from the circle due to the 
multiple (Coulomb) scattering. We calculate this function using the expression for the 
probability distribution function of muon lateral displacements based on the solution of 
a special kind of the adjoint transport equation obtained by [[|l[ . 

For solution of adjoint cascade equations we use a numerical method which is 



described in [42]. This method has been applied to solve a wide range of problems 



connected with cosmic ray physics (see for review |2"5|). The method [|2| provides an 
accuracy of solution of adjoint equations better than a few per cent and gives results 
for an arbitrary region (E Q , E mSuX ) of the primary energy in one run of calculations. 
The computational time consumed by the method is quite small and depends weakly 
(logarithmically) on the primary energy. For example, to calculate the number of cascade 
electrons (or photons) for a fixed altitude of the observation level and the primary energy 
range (10 9 , 10 22 ) eV one needs only a few minutes of the 1 GHz IBM PC computational 
time. 

The results for proton-induced air showers also have been obtained by solution of 
adjoint cascade equations with the help of the numerical method developed in Ref. 
|42j . Here we do not present the expressions for these equations, because they are quite 



similar to Eqs. (0) and (||). They can be found in fl23[ . 

4. Comparison with previous works 

A number of comparisons of air shower characteristics with similar results of previous 
studies has been performed to test the computational code. In figure 2 the cascade curves 
of photons are presented for the 7-ray induced showers developing in the homogeneous 
magnetic field (without presence of matter). One can see a perfect agreement with 
calculations of A.Goncharov 1271 . 



Figure 3 represents the cascade curves of electrons for air showers initiated by 
primary 7-rays of energy 10 20 eV, taking into account interactions with the geomagnetic 
field and the LPM effect. Three different combinations of the GMF and LPM effects 
are shown by curves 1 (no LPM, no GMF), 2 (only LPM), 3 (both LPM and GMF). A 
quite satisfactory agreement with the results of Ref. [^2] is seen for all 3 cases. 

In table 2 we present the number of muons at the sea level versus the total number 
of electrons for air showers created by primary 7-rays with energy below 10 16 eV 
(N e < 10 7 ). Our results are compared with similar results of Edwards et.al. [[TjJ. Note 
that the calculations of Ref. |19| , as well as many other studies in this energy region, 
were performed using simplified treatments of photohadronic interactions. Namely, it 
was assumed that the 7-hadron interaction is identical to the inelastic interaction of the 
pion of the same energy. Nevertheless, a quite satisfactory agreement is seen between 
results presented in the table. 



Table 2. The number of muons with energy T > 1 GeV as a function of the total 
number of electrons at the sea level for air showers created by inclined 7-rays with 
zenith angle 25° 



N e 


10 5 


3-10 5 


10 6 


3-10 6 


10 7 


Nfj,, present study 
N,*, @ 


0.29- 10 3 
0.33 -10 3 


0.77 -10 3 
0.92 • 10 3 


0.20 • 10 4 
0.27 -10 4 


0.51 • 10 4 
0.65 • 10 4 


0.15- 10 5 
0.16- 10 5 



In table 3 we present the calculation results for the number of photoproduced muons 
with kinetic energy T > 0. For comparison we show also the results of calculations 



performed by us using the CORSIKA code ||45f| . It should be noted that the results 
in table 3 correspond to a limited primary energy region which allows us to use the 
CORSIKA option based on the "complete" Monte Carlo method, i.e. without the 
application of the "thinning technique" f46fl . It is seen that discrepancy of results 
presented in the table does not exceed 20 %. 

The "complete" Monte Carlo method requires a very large computational time which 
does not allow extension of simulation to extremely high energy region. To reduce the 



computational time one has to apply the "thinning technique" [46]. Such calculations 



were performed by Capdevielle et al. [31] using the same CORSIKA simulation code for 



EHE primary 7-rays. In table 4 we compare our results on the muon content in EHE 



7-ray induced showers with similar results of Ref. [31 . Note that Capdevielle et al 



31] included in their calculations the LPM effect but ignored the effect of interactions 
of 7-rays with the geomagnetic field. Therefore in table 4 we show our results without 
the GMF effect as well. The results are in good agreement within 10 %. 



Table 3. The muon content in air showers created by primary 7-rays for z a ^ s =890 
g/cm 2 . Results of the numerical solution of adjoint equations are compared with 
simulation results obtained by the CORSIKA code. 

Primary energy of air shower, eV 3 • 10 13 10 14 3 • 10 14 10 14 3 • 10 14 

Zenith angle of arriving 7-ray, degree 45 45 

iV M , solution of adjoint equations 13.9 56.4 199 21.9 69.8 

N^, simulations with CORSIKA 11.6 50.2 174 26.1 79.0 



In figure 4 we present the radial distributions of the photoproduced muons for 
vertical air showers created by 7-rays with energy 10 15 eV. Results of the numerical 
solution of adjoint equations are compared with the simulations made by the "complete" 
Monte Carlo method using the CORSIKA code. [J. A perfect agreement of the results 
is achieved. 



Table 4. The number of muons with energy T > as a function of energy of 
the primary 7-ray for inclined (45°) air showers. The LPM effect is included, the 
geomagnetic field is disabled. 

EL, eV 10 17 10 18 10 19 10 20 



Nf,, present study 4.5 • 10 4 6.1 • 10 5 8.2 • 10 6 7.5 • 10 7 
Np, p} 5.0 -10 4 6.6 -10 5 8.1 -10 6 7.9 • 10 7 



In figure 5 we present the electron cascade curves for proton induced air showers. 
It is seen that our results are in agreement within 20 % with simulations performed in 
the framework of the QGS model [|4|. 



Table 5. The number of muons with energy T > for vertical air showers created by 
primary protons. Z b s =920 g/cm 2 . 



E p ,eV 10 18 3-10 18 10 19 3-10 19 10 20 

Nf, 7 present study 0.71 • 10 7 0.19 • 10 8 0.56 • 10 8 0.15 • 10 9 0.43 • 10 9 
N^, J47| 0.68 -10 7 0.20 -10 8 0.51 • 10 8 0.16 • 10 9 0.45 • 10 9 



In table 5 we present the muon content in air showers initiated by vertical protons. 
Our results on this quantity are compared with similar results of Nagano et.al. [47[ 



f The "thinning technique" |4q] does not allow adequate calculations of the radial distributions of 
muons in EHE electromagnetic air showers (see section 6), therefore the comparison presented in the 
table is limited by the primary energy 10 15 eV. 
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obtained by the CORSIKA simulation code with use of the QGSJET generator of hadron 
interactions. The difference between these two methods does not exceed 10 %. 



5. Longitudinal development of air showers 

5.1. Muon component of 7 -induced air showers 

In the energy region below ~ 10 16 eV the muon content in electromagnetic air showers 



has been studied by many authors (see e.g. |T8|, [19|, g(J, ^T], ^8], |49fl). However, the 
calculations , performed in this energy region, are based, as a rule, on rather simplified 
treatments of the photoproduction process. For example, it was often assumed that 
the 7-hadron interactions were identical to the inelastic pion interactions. Recently the 
photohadronic interactions have been implemented into the simulation codes CORSIKA 



[45|| and AIRES [|50| allowing accurate calculations of characteristics of muons in EHE 



electromagnetic air showers |5T> |36|, 51]. However, the interaction with the geomagnetic 



field has not yet been included in these codes. The latter reduces essentially the 
applicability of calculation results on the content of photoproduced muons calculated 
with the help of these codes. In the present study we calculate the shower characteristics 
using the technique of adjoint cascade equations and considering interactions of primary 
7-rays with a composite target consisting of two parts - the geomagnetic field and the 
Earth's atmosphere. 

In figure 6 we present the dependences of the total number of muons (-/V M ) at the 
observation level on the energy of primary 7-ray for vertical and inclined showers and 
for different assumptions about the LPM effect and the geomagnetic field. The following 
conclusions can be obtained from this figure. 

• If one uses the Bethe- Header's cross-sections for high energy interactions of cascade 
particles, the energy dependence of N^ is close to the power law, i.e. N^ ~ E& with 
/3~ 0.9 -1.0. 

• The LPM effect changes dramatically the energy dependence of N^. The growth 
of Nn(E^) becomes less rapid; moreover a local maximum arises near E 1 ~ 
10 20 — 10 21 eV, after which N^ decreases with E 1 (curves "2" in figures 6a and 
6b). 

• The GMF (without the LPM effect) does not essentially change the energy 
dependence of N^(E^). It remains to be continuously growing for both the vertical 
and inclined showers. 

• After including of the interactions with the geomagnetic field, the LPM effect does 
not affect considerably the A r At (£' 7 ) dependence. However, some deviation from 
the power law takes place in the energy region 10 19 -j- 10 21 eV. This deviation is 
connected with the fact, the LPM effects becomes effective a bit earlier (i.e. at 
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lower energies) compared with the GMF effect. 
• Behaviour of N^(E y ) for inclined air showers is generally similar to the one of the 
vertical showers. However, a local maximum of iV /i (£' 7 ) dependence, connected 
with the LPM effect, is observed for larger primary energies. Besides, for inclined 
showers the simultaneous influence of the LPM effect and the GMF can lead to a 
slight reduction of N^ in the high energy region, whereas for vertical air showers 
these effects lead to a growth of N^. 

The GMF influence on iV M as well as on other shower characteristics has a simple 
explanation. For sufficiently high product E n ■ B n the primary 7-ray cannot freely 
enter the Earth's atmosphere due to the pair production in the geomagnetic field. 
The secondary (pair-produced) electrons and positrons immediately interact with the 
same field through the synchrotron radiation. The free pathlength of these electrons is 
rather small (see figure 7) compared with the scale of the magnetosphere, and weakly 
depends on the particle energy. Thus, these electrons emit on their way to the Earth 
a large number of relatively low energy photons (according to figure 7 - one photon 
per every 10-100 km of the trajectory). On the other hand, the mean free pathlength 
for the pair production grows very rapidly with reduction of the photon energy (see 
figure 7). As a result, the energy of primary 7-ray is distributed among a large number 
of secondary photons with relatively small energies not able to interact effectively with 
the geomagnetic field (this effect is illustrated by table 6). Thus, instead of a single 
7-ray photon we have at the upper border of the atmosphere a 'bunch' of low energy 
photons (E < 10 19 eV) which thus appears below the effective threshold of the LPM 
effect. 



Table 6. The mean number of cascade synchrotron photons produced in GMF by a 
vertical primary 7-ray of energy 10 21 eV, and entering the Earth's atmosphere with 
energy greater than E t h- 

E th ,eV 10 17 10 18 10 19 



7V 7 430 135 28 



In figure 8 we present the integral energy spectra of photoproduced muons. These 
spectra correspond to different values of the energy and arrival direction of primary 
7-ray It is seen that for vertical air showers the spectrum is rather soft; only ~30% of 
muons have at the observation level kinetic energy exceeding 1 GeV. However, at large 
arrival angles the spectrum becomes considerably flatter. For example, for 9 = 45° 
about 40-45% of muons have energy larger than 1 GeV. At the same time, the energy 
spectrum of photoproduced muons rather weakly depends on the primary energy. 

In order to study the relation between the numbers of electrons and muons we 
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present in figure 9 the energy dependence of the ratio N^/N e for different arrival 
directions. It is seen that for vertical showers the total number of muons at the 
observation level is very small compared with the number of electrons. Both N e and 
N^ decrease with 6, but N e decreases faster. Thus at large zenith angles (6 > 65°) the 
number of muons becomes comparable or even can exceed the number of electrons. 

5.2. Comparison of muon content for air showers of different nature 

In figure 10 we present the energy dependence of the total number N^ of muons for 
vertical showers initiated by primary 7-rays and protons. The results of the present 



study are compared here with results for primary 7-rays obtained in Ref. [22] using 



an approximate analytical solution of cascade equations, and with the results of Ref. 



52 1 for a primary proton obtained by the Monte Carlo method. On the basis of 



comparison of data of Ref. j2l| and Ref. j|3) a conclusion has been made in Ref. [p2] 



that for EHE air showers the quantity N^ for 7-primary can be comparable with the 
number of muons in proton air showers. It is seen from figure 10 that the results of 
present calculations on A^ for 7-primary (curve 1) are by a factor of a few smaller than 
analytical results of []22| (curve 2). This discrepancy is explained by usage in Ref. [^ 



of a simplified model of the photoproduction process. It was assumed that the mean 
number of charged pions created in one 7-nuclear interaction did not depend on the 
7-ray energy (E^) and was equal to 2. Furthermore, the energy of each of these pions 
was assumed to be uniformly distributed in the interval (0,J5 7 ). This implies a very hard 
spectrum of secondary muons and an effective muon multiplicity of about 2. In fact, 
the multiplicity of the muons with the total energy more than 0.3 GeV [f] approaches 
to 2 only for parent 7-rays with energy equal to several GeV. It is demonstrated in 
figure 12 obtained using the SOPHIA generator ||38|| . On the other hand, the bulk of 
muons are produced in the electromagnetic air shower by cascade photons with typical 
energy of about 1 GeV for which the muon multiplicity is significantly less, n M ~ 0.5 
(see figure 12). At the same time, for determination of the ratio of muon contents in 7- 



and proton induced showers, Aharonian et al. |22[ used the calculations of McComb et 
al. [Q. As it is seen from figure 10, the calculations of ref. |43| (curve 5) in fact predict 
significantly smaller number of muons compared with our calculations (curve 3). Thus, 



the present study does not confirm the conclusion of Ref. |22| concerning the very high 
ratio of Nj?>/Nj?\ 

In figure 13 we present the calculation results on the energy dependence of the ratio 
5 = N^ /Nff) . One can see that this dependence has a non-trivial form. The behaviour 
of S can be explained by the fact the LPM effect starts "to operate" a bit earlier than 
the GMF does. In the energy region above 10 18 eV the ratio S ranges within ~ 0.1 -j- 0.2. 

f Exactly this value of the muon threshold energy is considered in figure 10. 
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In the energy region below 10 18 eV it decreases rapidly with reduction of the primary 
energy in accordance with the prediction of Ref. |22| . 



5.3. Electron component 

In figure 13 the dependence of the number of electrons on the energy of primary 7-ray 
is shown for the El Nuhuil site with use of different assumptions about the LPM effect 
and the geomagnetic field. One can see that the peculiarities of behaviour of N e are 
similar to the energy dependence of the muon content described in section 5.1 (compare 
figures 6 and 13). 

In figure 14 we show the results of calculations for the cascade curves of electrons. 
Such curves can be used for an approximate description of the temporal profile of the 
fluorescent radiation of air showers. In the case of 7-ray primary we show the cascade 
curves corresponding to both Northern and Southern Hemisphere sites of the Auger 
observatory. No significant difference between the shower profiles corresponding to 7- 
ray and proton induced showers can be seen. The difference in the GMF strength 
between Northern and Southern site locations (see table 1 and figure 1) does not affect 
considerably the shower profile as well. In this regard, an effective separation of 7-ray 
and proton-induced showers on the basis of the shower profiles seems to be hard. Perhaps 
such a separation could be achieved exploiting differences in the intrinsic fluctuations of 
the fluorescent light intensity, because the fluctuations in electromagnetic showers are 
smaller than in the hadronic ones (see e.g. |53|| ). The electromagnetic showers in the 
"photon bunch" regime caused by the GMF effect lead to an additional reduction of 
intrinsic fluctuations and, therefore, may improve the 7-ray separation efficiency 

6. Lateral distribution of photoproduced muons 

No detailed studies have been done until now on the lateral distribution function (LDF) 



of muons in electromagnetic EHE air showers. Recently, Capdevielle et al. |3T| published 
LDF of muons using the CORSIKA code and taking into account the LPM effect. 
However, they ignored the effect of interactions of 7-rays with GMF, which significantly 
reduces the applicability of these results. Moreover, the "thinning technique" used by 
these authors does not provide an acceptable accuracy of calculations for the LDF of 
muons. 



The thinning technique has been introduced by Hillas [46]. It allows a "complete" 
simulation only for a high energy part of the cascade which contains secondary particles 
with energy larger than a fixed threshold value E t ^. For low energy interactions, i.e. in 
the energy region < E t h, an artificial absorption of secondary particles is introduced. It is 
assumed that only one secondary particle survives in each interaction. For compensation 
of this " disappearance" of particles one should introduce appropriate statistical weights. 
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To avoid large fluctuations of these weights (dramatically worsening the accuracy of 
calculations), it is necessary to keep the thinning parameter e = E t h/E to be sufficiently 
small (E is the energy of the primary particle) . The computational time needed for 
simulation of one air shower by this technique is approximately proportional to e _1 . 



The simulations of [31] include the LDF of photoproduced muons for 10 individual 
7-ray induced air showers with primary energy 10 20 eV. These results correspond to the 
value of thinning parameter e = 10~ 6 . For such value of e one needs approximately 
one hour of the computational time of 500 MHz IBM PC to simulate one air shower 



45fl . It is seen from figure 7 of Ref. |3TJ that fluctuations caused by introduction of 
statistical weights are extremely large (up to a factor of 30). Therefore to provide an 
acceptably small statistical error (for example, 10%) one has to simulate an extremely 
large (up to ~ 10 5 ) number of air showers. Correspondingly, this requires a great deal 
of the computational time. For comparison, the calculation of LDF of photoproduced 
muons based on the numerical solution of adjoint cascade equations is not laborious and 
requires less than 1 hour of the computational time. Besides, this method gives a set of 
LDF in an arbitrary energy region of primaries just for one run of calculations. Note 
also that this method does not contain, by definition, statistical error at all. 

In figures 15, 16 we present the normalized LDF of muons of all energies, i.e. with 
kinetic energy T > T th = 0. We define this parameter as / M (r) = p^{r)/N^ where 
Pn(r) is the density of the muon flux at distance r from the shower axis. Results of 
figure 15 correspond to different assumptions about the GMF and the LPM effects, 
whereas figure 16 illustrates dependence of the LDF shape on the primary energy. It is 
seen that the geomagnetic field only slightly modifies the LDF shape. The interaction 
of 7-rays with the GMF makes the LDF shape insensitive to the LPM effect. Also, 
the LDF only slightly depends on the primary energy (see figure 16). For example, an 
increase of E 1 by a factor of 10 changes the LDF no more than ~20%. On the contrary, 
at the absence of GMF the LPM effect would steepen considerably the shape of LDF, 
and make it strongly dependent on the primary energy (figure 15). 

In figure 17 we present the muon density for vertical proton and 7-ray induced 
showers. We compare in this figure our results for primary 7-ray with the results for 



primary protons obtained in Ref. gy| by simulations with CORSIKA. One can see that 
the 7-induced shower provides an essentially narrower lateral distribution of muons 
compared with the proton shower. As a result, the ratio A = py / p$ decreases rapidly 
with a growth of the radius (from A ~ 0.1 at r = 200 m to A ~ 0.01 at r = 2500 m). 
This feature of A can provide an effective separation of 7-induced air showers by the 
Auger Observatory, for which the distance from the core of the shower to the most of 
triggered detectors is very large, several hundred meters or more. 

In figure 18 we illustrate dependence of muon LDF on the arrival direction of the 
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shower]]]. One can see that the width of the LDF increases rapidly with a growth 
of the zenith angle 9. This feature of LDF provides a significant contribution of the 
photoproduced muons to the total number of charged particles at large distances from 
the shower core. This effect of inclined showers is demonstrated in figure 19, where 
we present the particle density for two kinds of secondaries - muons and electronsf|. In 
our case the total number of muons at the observation level is much smaller than the 
total number of electrons (JV^ ~ 0.03iV e ; see figure 9). At the same time , at distances 
r > 1000 m from the core of the shower inclined at 65°, the muon density p M exceeds 
considerably the electron density p e . This effect may have a non- negligible impact on 
the estimates of the primary energy of inclined electromagnetic showers detected by 
air shower arrays not identifying the origin of secondary charged particles. At the 
same time, for the Auger Observatory allowing identification of secondary particles, the 
detection of unusually high /i/e ratio in inclined showers at large distances from the 
shower core should not be a priori accepted as hadronic showers. 

7. Conclusion 

The technique of adjoint cascade equations has been applied to study the characteristics 
of air showers created by EHE 7-rays such as the total numbers of cascade electrons and 
photoproduced muons at the observation level, the longitudinal and lateral distributions 
of electrons and muons. Both the LPM effect and interactions of primary 7-rays with 
the geomagnetic field have been incorporated in the computational code. 

The LPM effect changes dramatically the energy dependence of the muon content 
Nfj,. Due to this effect the growth of N^ with energy gradually disappears producing a 
local maximum around _E 7 ~ 10 20 — 10 21 eV, after which iV M decreases with E y . The 
interactions of 7-rays with GMF which start to "operate" almost simultaneously with 
the LPM effect (i.e. at E > 10 19 eV) recover the nominal energy dependence of iV M 
"switching off" the LPM effect. 

The ratio of muon contents in the 7-ray and proton induced showers ranges between 
0.1 — 0.2 in the EHE energy region, but decreases significantly at lower energies. The 
energy dependence of this quantity has a rather complicated (non power-law) form. 

Due to the combination of LPM and GMF effects the cascade curves of electrons in 
7-ray induced showers are quite similar to these curves for the proton-induced showers. 
This makes unfortunately rather difficult the separation of 7-induced air showers on the 
basis of the temporal profiles of the fluorescent light by the Fly's Eye type detectors, 
except a narrow primary energy band (around 10 19 eV) and a specific region of arrival 

f For the inclined air showers we consider the LDF averaged over the azimuth angle in a plane 

perpendicular to the shower arrival direction. 

J For the LDF of electrons we use data from Ref. [H . 
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directions (depending on the detector location site) for which 7-rays freely pass through 
the Earth's magnetosphere. 

The effect of LPM, when combined with interactions of 7-rays with the GMF, on the 
lateral distribution of photoproduced muons is quite small. The 7-induced air shower 
has an essentially narrower lateral distribution of muons compared with the proton 
shower. This feature can be effectively used for separation of 7-ray induced air showers 
by the particle detectors of the Auger Observatory. The extension of the muon LDF 
increases considerably with the zenith angle of arriving primary 7-rays. This results 
in an interesting effect - the density of muons at large distances from the shower core 
(> 1000 m) may noticeably exceed the electron density. 
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Figure 1. Profiles of the perpendicular component of geomagnetic field for 
different detector locations and for different zenith angles of arriving 7-ray 
photon. 
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Figure 2. The number of secondary synchrotron photons with energies 
E > 10 17 , 10 18 and 10 19 eV produced in the homogeneous magnetic field 
with B n — 0.3 G by primary 7-rays of energy 10 21 eV. Curves - this work; 
points - from Ref. |27l . 
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Figure 3. Dependence of the number of cascade electrons with energy 
T > on the depth of the observation level for vertical primary 7-rays of 
energy 10 20 eV with use of different assumptions concerning the LPM and 
GMF effects. Solid curves - this work; dashed curves - from Ref. |E2|. 
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Figure 4. The radial dependence of the density of muons with energy T > 
for vertical air showers produced by primary 7-rays with energy E = 10 15 eV 
at the sea level (1) and at 5 km above the sea level (2). Curves - numerical 
solution of adjoint equations; points - simulations with CORSIKA. 
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Figure 5. The number of electrons with energy T > in the air showers 
initiated by vertical protons as a function of the penetration depth. The 
energies of primary protons are indicated at the curves. Solid curves - 
present work; dashed curves - from Ref . H] . 
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Figure 6. Dependence of the number of muons with energy T > on the 
energy of primary 7-ray for different assumptions concerning the LPM and 
GMF effects calculated for the El Nuhuil site. Vertical (a) and inclined ( b, 
45°) showers are considered. 
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Figure 7. The energy dependence of the mean free pathlength of 7-rays 
and electrons in the homogeneous magnetic field B n = 0.3 and 0.03 G. 
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Figure 8. The fraction of muons with kinetic energy larger than T t h for 
different incident angles 9 and energies E 1 of primary 7-ray. The El Nuhuil 
site. 
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Figure 9. The ratio of total numbers of muons and electrons as a function 
of energy of primary 7-ray. The zenith angles of the arriving 7-rays are 
indicated at the curves. The El Nuhuil site. 
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Figure 10. The energy dependence of number of muons with total energy 
E t h > 0.3 GeV for vertical 7-ray and proton induced showers at the sea 
level. 
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Figure 11. The mean multiplicity of muons with kinetic energy T > T t h 
produced in the 7-nuclcar interactions as a function of the 7-ray energy. The 
values of T t h are indicated at the figures. 
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Figure 12. The ratio of the contents of muons with energy T > for 
air showers produced by primary 7-rays and protons. The arrival angles of 
primary particles are indicated at the curves. The El Nuhuil site. 
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Figure 13. Dependence of the number of cascade electrons with energy 
T > on the energy of vertical 7-ray and proton induced air showers for 
different assumptions concerning the LPM and GMF effects. The El Nuhuil 

site. 
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Figure 14. The number of cascade electrons with energy T > as a 
function of the penetration depth for primary protons and 7-rays of energy 
10 20 eV. Two different detector locations - El Nuhuil and Millard County - 
are considered for inclined (60°) air showers. 
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Figure 15. The lateral distributions of muons with energy T > for 
vertical showers induced by a primary 7-ray of energy 10 20 eV, calculated 
for different assumptions concerning the LPM and GMF effects. The El 
Nuhuil site. 
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Figure 16. The lateral distributions of muons with energy T > for vertical 
7-ray showers of energies indicated at the curves. The El Nuhuil site. 
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Figure 17. The muon density in the vertical showers produced by protons 
and 7-rays of energy 10 20 eV. The results of present work for primary 7- 
ray (the El Nuhuil site) and for three different muon threshold energies 
are compared with results of Ref. 1 )4 7| calculated for primary protons at 
Z oha = 920 g/cm 2 . 
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Figure 18. The lateral distributions of muons with energy T > in the 
electromagnetic showers for 3 different arrival angles of primary 7-ray of 
energy 10 20 eV. The El Nuhuil site. 
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Figure 19. The densities of muons (p M ) and electrons (p e ) in the inclined 
(65°) 7-ray shower of energy 10 20 eV. The El Nuhuil site. 



